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Abstract 

The  purpose  of  this  research  is  to  establish  structure-forming  principles  that  govern  the  atomic  structures  of  metallic  glasses.  A structural 
model  based  on  efficient  atomic  packing  will  be  discussed  and  applied  to  the  topological  systems  that  represent  most  metallic  glass  alloys. 
The  concept  of  efficient  atomic  packing  has  direct  and  specific  implications  regarding  the  local  structure  and  composition  of  metallic  glasses. 
Specific  solute-to-solvent  atomic  radius  ratios  and  specific  solute  concentrations  related  to  these  ratios  are  shown  to  be  preferred  in  this  model, 
and  analysis  of  a wide  range  of  metallic  glass  systems  shows  a very  strong  correlation  with  these  predicted  values.  Relationships  between 
atomic  size  and  concentration  are  discussed,  and  new  insights  are  proposed  based  on  the  current  structural  model.  Possible  local  atomic 
configurations  (i.e.,  atomic  clusters)  are  defined  based  on  topological  constraints  that  are  derived  from  the  requirement  of  efficient  atomic 
packing.  Experimental  observations  drawn  from  the  literature  that  provide  support  for  this  model  are  presented. 

© 2003  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

A number  of  structural  theories  have  been  developed 
for  the  formation  of  glasses  [1,2].  The  most  widely  known 
theory  is  attributed  to  Zachariason  [3]  and  provides  four 
rules  that  describe  the  formation  and  linking  of  character- 
istic coordination  polyhedra.  This  model  offers  a robust 
framework  for  describing  the  structure  of  oxide  glasses, 
including  those  based  on  silicates,  borates  and  phosphates. 
Bonding  in  these  materials  includes  significant  amounts 
of  ionic  (non-directional)  and  covalent  (directional)  bonds. 
Less-known  structural  theories  include  Sun’s  single  bond 
strength  criterion  [4],  based  on  bond  strength,  and  Dietzel’s 
field  strength  criterion  [5],  based  on  the  magnitude  of 
electrostatic  interaction.  More  recently,  a structural  theory 
has  been  developed  for  chalcogenide  glasses  [6].  Based 
on  topological  considerations,  this  model  provides  insights 
into  glasses  formed  in  solids  that  are  dominated  by  covalent 
bonding.  The  structures  formed  in  all  of  these  glass  systems 
are  relatively  inefficiently  packed  due  to  constraints  aris- 
ing from  preferred  bond  angles  in  solids  with  a significant 
covalent  character  and  from  constraints  that  arise  from  the 
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requirement  to  maintain  charge  neutrality  in  ionic  solids. 
Thus,  the  coordination  number  ( N)  is  low  in  these  glass 
systems  and  is  generally  less  than  ~8. 

In  contrast,  cohesion  in  metallic  glasses  is  dominated  by 
metallic  bonding.  Although  covalent  bonding  is  sometimes 
suggested  for  metal-metalloid  [7]  and  Al-transition  metal 
(TM)  [8,9]  glasses,  it  has  been  argued  that  strong  covalent 
bonding  is  difficult  to  rationalize  for  the  metal-metalloid 
glasses  due  to  the  large  coordination  numbers  (8-9)  and  the 
absence  of  boron  d-orbitals  [1].  By  removing  constraints  in- 
troduced by  preferred  bond  angles  in  covalent  solids  and,  to 
a lesser  extent,  local  charge  neutrality  in  ionic  solids,  atomic 
packing  in  metallic  glasses  can  be  very  efficient.  Thus, 
atomic  packing  in  crystalline  metals  is  generally  higher  than 
in  solids  with  a significant  degree  of  covalent  bonding,  such 
as  silicon  and  diamond.  Furthermore,  the  density  of  the  best 
metallic  glasses  is  an  exceptionally  high  fraction  (~0.995) 
of  the  density  of  the  same  alloy  in  the  crystalline  state 
[10,1 1],  As  a result,  N in  metallic  glasses  is  rarely  less  than 
8,  and  can  be  as  large  as  17  [12].  In  spite  of  theories  based 
on  randomness  as  a structural  paradigm  for  metallic  glasses, 
significant  short-range  chemical  and  topological  ordering  is 
common  in  metallic  glasses  (see,  for  example,  [9,13,14]). 
The  general  features  described  here  for  the  structure  of 
metallic  glasses  persist  even  in  glasses  with  a small  amount 
of  covalent  bonding  in  addition  to  metallic  bonding.  Thus, 
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a structural  theory  for  the  formation  of  metallic  glasses  that 
is  able  to  explain  the  important  structural  distinctions  of 
coordination  number  and  atomic  packing  efficiency  relative 
to  oxide  and  chalcogenide  glasses  is  required. 

Earlier  efforts  to  develop  a structural  theory  for  metallic 
glasses  have  considered  concepts  built  upon  a dense  ran- 
dom packing  (DRP)  paradigm.  Efforts  to  rationalize  the  high 
relative  density  of  glasses  in  the  amorphous  structure  rela- 
tive to  the  same  alloy  composition  in  the  crystalline  state 
(prel  = Pamor/Pxtai)  have  emphasized  placement  of  solute 
atoms  in  the  inefficiently  packed  local  regions  (‘holes’)  that 
characterize  a DRP  array  of  spheres  [15-17],  Careful  analy- 
ses have  shown  that  these  approaches  do  not  provide  quan- 
titative agreement  with  the  high  prei  of  metallic  glasses 
[17-19],  A different  approach,  based  on  a continuous  linked 
array  of  randomly  arranged  trigonal  prismatic  clusters  in 
metal-metalloid  glasses,  has  been  proposed  but  has  not  been 
fully  developed  [13]. 

The  purpose  of  this  manuscript  is  to  summarize  and  dis- 
cuss recent  developments  in  the  establishment  of  structure- 
forming principles  that  govern  the  atomic  structure  of 
metallic  glasses.  Efforts  to  establish  a structural  paradigm 
based  on  efficient  atomic  packing  will  be  described,  and  the 
influence  of  system  topology  (relative  sizes  and  concentra- 
tions of  constituent  atoms)  on  local  atomic  configurations 
will  be  discussed.  Initial  models  have  emphasized  the  rela- 
tive sizes  of  constituent  atoms.  The  primary  metrics  used  to 
determine  the  validity  of  these  models  will  include  consis- 
tency with  widely  accepted  physical  principles,  an  ability 
to  explain  the  basic  features  for  metallic  glasses  described 
above,  and  agreement  with  relevant  experimental  obser- 
vations related  to  the  formation  of  metallic  glasses.  It  is 
anticipated  that  establishment  of  such  principles  will  pro- 
vide an  improved  understanding  of  the  stability  of  metallic 
glasses,  and  may  lead  to  a new  structural  theory  for  the 
formation  and  stability  of  metallic  glasses. 

2.  The  principle  of  efficient  atomic  packing 

As  described  above,  the  best  metallic  glasses  possess  a 
high  density  relative  to  the  crystalline  state  of  the  same 
alloy,  which  suggests  that  efficiently  packed  atomic  con- 
figurations must  be  common  features  in  the  amorphous 
structure.  From  a kinetic  perspective,  an  efficiently  packed 
atomic  structure  is  expected  to  provide  higher  viscosity 
relative  to  a more  open  structure.  This  increased  viscos- 
ity reduces  the  rate  of  mass  transport,  and  hence  restricts 
nucleation  and  growth  of  competing  crystalline  phases. 
From  an  energetic  perspective,  a decrease  in  system  volume 
decreases  system  energy,  thus  stabilizing  the  amorphous 
structure  relative  to  competing  structures.  Consistent  with 
this  concept,  evidence  shows  that  metallic  glass  alloys  of- 
ten display  a minimum  in  molar  volume  at  the  best  glass 
forming  composition  [20,21].  Thus,  the  general  concept 
of  efficient  atomic  packing  has  a sound  physical  basis  that 


is  consistent  with  both  energetic  and  kinetic  models  of 
metallic  glass  formation.  While  a structural  model  based  on 
efficient  atomic  packing  must  consider  a range  in  relevant 
length  scales,  the  present  developments  will  emphasize  lo- 
cal packing  efficiency  on  a length  scale  of  a few  atomic 
diameters.  In  the  following  section,  the  influence  of  relative 
atomic  size  on  local  packing  efficiency  will  be  described. 

3.  Influence  of  relative  atomic  size 

A simple  model  is  presented  here  to  establish  the  influ- 
ence of  relative  atomic  size  on  local  packing  efficiency  and 
to  quantify  local  packing  efficiency.  This  will  be  described 
briefly  in  two  dimensions  (2D),  and  will  then  be  extended  to 
three  dimensions  (3D).  Consider  a single  circle  surrounded 
by  N equally  sized  larger  circles,  with  a radius  ratio  R be- 
tween the  central  (solute)  and  outer  (solvent)  circles.  If  R is 
selected  so  that  there  are  no  gaps  between  the  circles  in  the 
first  coordination  shell,  then  it  can  be  seen  that  this  is  qual- 
itatively an  efficiently  packed  configuration.  If  the  central 
circle  is  enlarged  slightly,  then  R increases  afid  N remains 
unchanged,  but  gaps  are  formed  between  the  solvent  circles, 
forming  a less  efficiently  packed  configuration.  As  R con- 
tinues to  increase,  the  gaps  between  solvent  circles  continue 
to  increase,  producing  a continuing  decrease  in  packing 
efficiency  ( P ).  Eventually,  a value  of  R is  reached  where  an 
efficiently  packed  configuration  with  N'  ~ N + 1 is  once 
again  achieved.  A simple  relationship  between  R and  the 
theoretical  coordination  number  ( Nr ) is  given  as  [22] 


arcsin(l/(l  + /?)) 

NJ  is  a real  number  where  the  integer  portion  represents  the 
number  of  full  coordinating  circles,  N,  and  the  fractional 
portion  is  a quantitative  measure  of  the  gaps  between  the 
circles  in  the  first  coordination  shell.  P may  then  be  defined 
as  the  actual  value  of  N normalized  by  NT.  It  can  be  seen 
from  this  discussion  that  the  local  packing  efficiency  in  2D 
will  reach  a maximum  value  of  unity  at  those  specific  values 
of  R where  ,¥T  is  an  integer. 

A similar  analysis  has  recently  been  extended  to  3D 
by  deriving  a solution  for  NJ  that  describes  the  number 
of  equal-sized  (solvent)  spheres  that  can  fit  around  a cen- 
tral (solute)  sphere  [19].  This  represents  a solute-centered 
atomic  cluster  with  solvent  atoms  only  in  the  first  coordi- 
nation shell.  The  result  is  given  as  [19] 


n(2  — g)+(2g,)arccos{(sin^/g)[l  - 1 /(R  + l)2]1/2} 

(2) 

where  q is  the  maximum  number  of  solvent  atoms  in  the  first 
coordination  shell  that  can  contact  a given  solvent  atom  that 
is  also  in  the  first  coordination  shell.  The  q value  depends 
explicitly  on  R and  is  3, 4 or  5.  A may  be  determined  exper- 
imentally from  diffraction  data,  or  may  be  calculated  as  the 
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Table  1 

Values  of  N and  corresponding  values  of  R* 


N 

R*  (Eq.  (2)) 

R*  (cluster) 

N 

**  (Eq.  (2)) 

R*  (cluster) 

3 

0.1 55a 

0.155 

12 

0.902 

0.902 

4 

0.225 

0.225 

13 

0.976 

ND 

5 

0.362 

Not  stable 

14 

1.047 

ND 

6 

0.414 

0.414 

15 

1.116 

ND 

7 

0.518 

0.591 

16 

1.183 

ND 

8 

0.617 

0.645 

17 

1.248 

ND 

9 

0.710 

0.732 

18 

1.311 

ND 

10 

0.799 

0.834 

19 

1.373 

ND 

11 

0.884 

0.902 

20 

1.433 

ND 

a Calculated  from  Eq.  (1).  ND:  not  determined. 


truncated  value  of  N7.  In  a hard  sphere  cluster,  such  as  those 
described  in  the  following  section,  N is  the  number  of  full 
spheres  in  the  first  coordination  shell.  P may  be  determined 
in  3D  by  truncating  N[  in  Eq.  (2)  and  dividing  by  Nr . As 
described  for  the  2D  example,  a maximum  P is  obtained  for 
specific  values  of  R where  NT  is  an  integer,  designated  as 
R*.  These  values  are  given  in  Table  1 for  a range  in  R that  is 
most  relevant  for  metallic  and  oxide  glasses.  Additional  de- 
tails relating  to  the  derivation  of  these  values  is  given  in  [19], 

Thus,  if  efficient  atomic  packing  in  the  first  coordination 
shell  is  important  in  the  formation  of  metallic  glasses,  then 
glasses  should  show  a preference  for  these  solute-to-solvent 
atom  radius  ratios.  A large  number  of  binary  and  complex 
metallic  glass  systems  has  been  analysed  [19]  to  explore  a 
possible  correlation  between  actual  atomic  radius  ratios  in 
reported  metallic  glasses  and  the  predicted  values  of  R*  in 
Table  1.  One  radius  ratio  (R  = rsoiute/rSolvent)>  representing 
solute-centered  clusters  with  solvent  atoms  only  in  the  first 
coordination  shell,  was  determined  for  each  of  76  binary 
metallic  glasses.  Complex  metallic  glass  alloys  based  on  Al, 
Au,  Fe,  La,  Mg,  Nd,  Ni,  Pd,  Pt,  Si,  Sm,  Ti  and  Zr  were  also 
included  in  this  analysis.  One  radius  ratio  was  determined  for 
each  solute  in  a given  alloy  family.  For  example,  five  values 
of  R were  determined  from  the  Zr-Ti-Al-Cu-Ni-Be  metal- 
lic glass  system — one  R for  each  of  the  five  solutes.  While 
atomic  radii  depend  upon  local  structure  and  chemistry,  dif- 
ferences are  generally  small  in  systems  where  metallic  bond- 
ing dominates.  Important  differences  from  handbook  values 
for  metallic  radii  exist  in  some  systems,  including  Al-TM 
glasses  and  glasses  containing  metalloids  such  as  B,  C,  or 
P.  The  atomic  radii  used  in  the  present  analysis  were  taken 
from  a critical  assessment  that  accounts  for  these  differences 
[23].  Additional  details  of  this  analysis  of  the  correlation 
between  solute  radius  ratios  and  predicted  values  of  R*  is 
provided  in  [19]. 

A summary  histogram  showing  the  frequency  with  which 
particular  radius  ratios  occur  within  these  metallic  glasses  is 
shown  in  Fig.  1 [19].  The  vertical  lines  represent  the  values 
of  R*  calculated  from  Eq.  (2).  A strong  correlation  is  shown 
for  each  class  of  metallic  glass,  and  for  all  metallic  glasses 
taken  together.  Although  not  yet  rigorously  analysed,  initial 
observations  of  solute-to-solvent  radius  ratios  in  represen- 


tative silicate  and  borosilicate  oxide  glasses  also  show  a 
correlation  with  the  predicted  values  of  R*  [19],  as  do  heavy 
metal  fluoride  glasses.  Based  on  the  strength  of  these  cor- 
relations, it  is  concluded  that  the  model  for  efficient  atomic 
packing  in  the  first  coordination  shell  is  a reasonable  phys- 
ical representation  of  the  local  atomic  structure  in  metallic 
glasses. 

A primary  feature  of  this  model  is  the  presence  of  solvent 
atoms  only  in  the  first  coordination  shell.  This  is  an  accurate 
representation  of  the  structure  of  silicate  and  borosilicate 
glasses.  In  addition,  there  is  essentially  no  metalloid- 
metalloid  bonding  in  metal-metalloid  glasses  [14,24],  so 
that  this  model  also  fairly  represents  these  well-known 
metallic  glass  systems.  While  significant  short-range  chem- 
ical ordering  is  a common  feature  in  many  metallic  glasses, 
some  degree  of  solute-solute  bonding  will  be  expected 
in  metallic  glasses.  Where  data  is  available,  in  fact,  some 
solute-solute  bonding  has  been  observed  [12,25-29].  Fur- 
ther, the  present  model  provides  a hard  sphere  idealization 
with  no  distortion  in  solute-solvent  bonds.  In  spite  of  these 
details,  the  strong  correlation  between  predicted  radius  ra- 
tios and  those  commonly  observed  in  a diverse  range  of 
metallic  glass  systems  suggests  that  these  deviations  in  real 
systems  from  the  model  idealization  may  represent  only  a 
small  correction  to  the  model  summarized  here. 

3.1.  Characteristic  atomic  configurations 

The  model  summarized  in  the  preceding  section  is  based 
on  the  maximum  number  of  solvent  atoms  that  may  exist 
in  the  first  coordination  shell  of  a given  solute  atom.  This 
is  accomplished  by  determining  the  minimum  area  on  a 
solute  atom  surface  that  is  associated  with  (or  ‘occupied 
by’)  a solvent  atom  in  the  first  coordination  shell.  The  total 
surface  area  of  the  solute  atom  is  divided  by  this  minimum 
area  associated  with  a single  solvent  atom  to  determine  NT. 
However,  the  actual  area  associated  with  a solvent  atom  is 
sometimes  larger  than  this  minimum  value.  For  example, 
an  octahedron  exhibits  N — 6 for  R = 0.414,  and  only 
equilateral  triangular  faces  are  produced  by  connecting  all 
of  the  vertices.  This  represents  the  most  efficient  packing 
around  the  central  solute  atom.  On  the  other  hand,  a trig- 
onal prism  also  possesses  N = 6,  but  R = 0.528  and  the 
faces  of  this  polyhedron  consist  of  two  equilateral  triangles 
and  three  squares.  The  three  square  faces  indicate  a less  ef- 
ficient packing,  so  that  the  surface  area  of  the  central  atom 
associated  with  each  of  the  atoms  in  the  first  coordination 
shell  is  larger  than  for  the  octahedron.  This  larger  surface 
area  associated  with  each  solvent  atom  results  in  a larger 
R for  the  trigonal  prism  configuration  relative  to  the  octa- 
hedron. Thus,  the  actual  value  of  R required  to  support  a 
given  value  of  N also  depends  upon  geometrical  constraints 
arising  from  packing  of  the  solvent  spheres  in  the  first  coor- 
dination shell.  As  a result,  the  actual  atomic  configurations 
must  be  specified  to  accurately  determine  the  relationship 
between  R and  N.  Stated  differently,  a given  value  of  N 
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(a)  Radius  Ratio,  R (b)  Radius  Ratio,  R 


Fig.  1.  Histograms  displaying  the  correlation  between  predicted  radius  ratios  (/?*),  shown  as  vertical  lines  in  the  figures,  and  solute-to-solvent  radius 
ratios  observed  in  (a)  binary  metallic  glasses,  (b)  metallic  glasses  with  marginal  glass  forming  ability  and  (c)  bulk  metallic  glasses.  A summation  of 
these  three  plots  is  shown  in  (d)  [19]. 


may  be  produced  by  several  values  of  R,  depending  on  the 
actual  configuration  of  atoms  in  the  first  coordination  shell. 
Alternatively,  NJ  represents  an  idealized  configuration,  so 
that  a unique  value  of  NJ  is  obtained  for  a given  value  of  R. 

Consideration  of  local  atomic  configurations  is  generally 
based  on  polyhedra  that  are  common  in  crystalline  or  qua- 
sicrystalline solids,  with  N — 4,  6,  8 or  12.  In  addition,  a 
capped  trigonal  prism  with  N = 9 has  been  established  as 
a representative  cluster  in  metal-metalloid  glasses  [13].  In 
metallic  glasses,  values  of  R are  rarely  smaller  than  ~0.6,  so 
that  N > 8 (see  Table  1).  Thus,  atomic  clusters  with  N = 10 
or  1 1 may  exist  in  metallic  glasses  with  appropriate  values 
of  R , but  the  atomic  configurations  possible  for  these  values 
of  N are  not  yet  established.  Further,  inspection  of  clusters 
with  N = 8 from  crystalline  solids  suggests  intuitively  that 
other  more  efficiently  packed  configurations  may  be  possi- 
ble. Relatively  little  work  has  been  reported  for  clusters  that 
are  likely  in  systems  with  R > 1.  In  Frank-Kasper  phases, 
atomic  clusters  with  N = 14,  15  or  16  have  been  described 
in  detail  and  have  relevance  in  topologically  close  packed 
structures  [30,31].  However,  these  clusters  were  constructed 
with  constraints  that  may  not  represent  the  most  efficiently 


packed  configurations.  Further,  solutes  with  R as  large  as  1.4 
are  often  observed  in  metallic  glasses,  so  that  clusters  with 
N as  large  as  19  may  be  relevant  in  metallic  glasses. 

A systematic  effort  to  explore  and  establish  the  actual 
atomic  configurations  that  may  exist  for  8 < N < 19  is 
now  underway,  and  the  initial  observations  are  presented 
here.  More  than  one  configuration  exists  for  most  values 
of  N,  and  these  different  configurations  result  in  a range 
of  R values  for  a given  N.  For  example,  six  distinct  cluster 
configurations  have  been  constructed  with  N = 10,  and 
corresponding  values  of  R range  from  0.834  to  0.902.  Only 
stable  configurations  are  considered,  defined  here  as  one 
where  all  equal-sized  spheres  in  the  first  coordination  shell 
just  touch  the  central  sphere,  and  each  sphere  in  the  first 
coordination  shell  touches  at  least  three  other  spheres  in  the 
first  shell.  After  constructing  possible  clusters,  R can  be  cal- 
culated for  each  specific  configuration  from  geometric  anal- 
ysis. The  minimum  value  of  R for  each  value  of  N provides 
a cluster-based  determination  of  R *,  representing  the  most 
efficiently  packed  configuration  for  that  value  of  N.  The 
cluster-based  values  for  R*  determined  in  this  way  are  shown 
in  Table  1,  along  with  the  values  calculated  from  Eq.  (2). 
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Fig.  2.  Candidate  clusters  with  high  packing  efficiency:  (a)  a distorted  capped  trigonal  prism  with  N = 9 and  R = 0.732,  and  (b)  a cluster  with  N = 10 
and  R = 0.834. 


Cluster-based  values  of  R*  for  N = 3,  4,  6 and  12  match 
the  idealized  predictions  of  Eq.  (2)  exactly.  With  the  excep- 
tion of  the  cluster  for  N = 1,  the  remaining  actual  clusters 
with  N < 12  have  values  of  R*  that  are  less  than  5%  higher 
than  the  ideal  values  predicted  from  Eq.  (2).  Studies  in  the 
mathematics  community  have  provided  rigorous  solutions 
for  the  minimum  value  of  R required  for  2 < N < 12  (see 
Table  11.1  in  [32]).  The  results  of  these  efforts  match  the 
cluster-based  values  of  R*  in  Table  1 to  four  significant  dig- 
its, validating  this  cluster-based  approach  for  determining 
values  of/C.  The  most  efficiently  packed  cluster  for  /V  = 8 
is  a square  antiprism  and  for  N = 1 1 is  an  icosahedron  with 
one  sphere  removed  from  the  first  coordination  shell.  Effi- 
ciently packed  clusters  with  N — 9 and  10  are  shown  in 
Fig.  2.  The  cluster  with  N — 9 is  a capped  trigonal  prism 
that  has  been  extended  along  the  prism  axis  to  allow  the 
three  ‘capping’  atoms  to  contact  the  solute  sphere  as  first 
nearest  neighbors  [33],  This  configuration  is  stable  for  R = 
0.732,  whereas  the  trigonal  prism  is  stable  for  R = 0.528. 
These  larger  values  of  R allow  the  distorted  capped  trig- 
onal prism  to  dominate  in  metal-metalloid  glasses,  where 
the  solute-to-solvent  radius  ratios  range  from  0.609  to  0.781 
[33], 

A systematic  investigation  of  configurations  that  will 
produce  a minimum  value  of  R for  13  < IV  < 19  is  also 
suggested  to  be  of  interest  with  respect  to  possible  atomic 
configurations  in  metallic  glass  structures.  Initial  work  sug- 
gests that  the  differences  in  cluster-based  values  of  R*  and 
those  from  Eq.  (2)  may  be  more  significant  than  for  N > 12. 
While  there  may  be  some  interest  in  simply  defining  such 
configurations  from  a fundamental  perspective,  it  is  also  im- 
portant to  identify  approaches  for  determining  whether  such 
configurations  actually  occur  in  metallic  glasses.  Critical 
observations  in  diffraction  experiments  may  provide  some 
evidence,  and  careful  analysis  of  atomistic  simulations  may 
also  provide  insight  and  support  for  such  configurations. 


Possible  idealized  local  atomic  configurations  may  be  de- 
fined in  a particular  alloy  system  by  considering  efficiently 
packed  clusters  that  may  occur  at  the  particular  radius  ra- 
tios exhibited  by  that  alloy  system.  Support  for  the  presence 
of  such  local  configurations  may  be  explored  by  compar- 
ing model  predictions  based  on  coordination  numbers  and 
measurement  of  local  packing  efficiency  with  experimental 
data.  This  approach  has  recently  been  pursued  in  Al-Y  and 
Al-Y-Ni  glass  alloys  [34].  The  results  of  this  initial  study 
are  consistent  with  the  experimental  data,  so  that  the  basic 
features  of  this  model  may  be  relevant  in  the  actual  metallic 
glass.  However,  further  investigation  is  required  to  provide 
more  convincing  evidence  for  the  validity  of  this  model.  Ad- 
ditional investigations  may  include  comparisons  with  higher 
quality  diffraction  data  than  those  currently  available  in  the 
literature  and  may  also  include  analysis  of  atomic  simula- 
tions. Such  efforts  are  now  underway. 


4.  Influence  of  relative  atomic  concentration 

Earlier  efforts  to  establish  the  importance  of  atomic  con- 
centrations of  constituent  atoms  on  glass  stability  have  fo- 
cused on  the  relationship  between  atomic  size  and  atomic 
concentration.  A well-known  model  [23]  provides  an  explicit 
relationship  relating  the  volume  difference  between  solute 
and  solvent  atoms  and  the  minimum  solute  concentration  re- 
quired to  destabilize  the  competing  binary  crystalline  phase. 
This  model  gives  an  inverse  relationship  between  R and  the 
critical  solute  concentration  in  binary  glasses,  and  provides  a 
remarkably  good  ability  to  predict  the  critical  solute  concen- 
tration for  a significant  number  of  binary  metallic  glasses. 

A phenomenological  approach  for  visualizing  and  com- 
paring the  topology  of  binary  and  complex  metallic  glasses 
has  recently  been  devised  by  plotting  the  atomic  concen- 
tration and  the  atomic  size  relative  to  the  solvent  atom  for 


D.B.  Miracle  et  al/ Materials  Science  and  Engineering  A 375-377  (2004)  150-156 


155 


each  element  in  a metallic  glass  [35].  These  atomic  size  dis- 
tribution plots  (ASDPs)  have  been  constructed  for  a large 
number  of  binary  [33]  and  complex  [35,36]  metallic  glasses. 
While  ASDPs  for  most  metallic  glasses  produce  a profile 
of  R versus  concentration  that  is  consistent  with  Egami’s 
model,  most  of  the  best  metallic  glasses  show  a relationship 
where  solute  concentration  increases  as  R decreases  below  a 
value  of  0.8-0.9.  Similarly,  the  most  concentrated  solute  in 
silicate  glasses  (Si)  is  generally  the  smallest  solute.  Egami’s 
model  is  not  expected  to  predict  the  detailed  profiles  rep- 
resented in  ASDPs,  since  the  former  is  developed  for  bi- 
nary glasses  only  while  ASDPs  of  BMGs  represent  complex 
metallic  glasses  where  each  of  the  constituents  have  been 
optimised  to  improve  glass  formability.  Nevertheless,  the  in- 
ability of  Egami’s  model  to  represent  the  general  trend  in 
the  relationship  between  concentration  and  atomic  size  sug- 
gests that  some  fundamental  physical  feature  may  be  miss- 
ing from  this  model. 

While  Egami’s  model  considers  only  substitutional  occu- 
pancy for  solute  atoms,  recent  modifications  consider  that 
solute  atoms  may  occupy  either  substitutional  or  interstitial 
sites  in  the  competing  crystalline  lattice  [37,38].  For  R < 
0.81,  interstitial  occupancy  is  energetically  favored  (pro- 
duces a smaller  strain  energy)  over  placing  the  same  solute 
in  a substitutional  site  [38],  Further,  as  R becomes  increas- 
ingly smaller  than  0.81,  the  strain  associated  with  an  intersti- 
tial solute  becomes  increasingly  smaller,  so  that  the  critical 
concentration  increases  with  increasing  difference  between 
solute  and  solvent  atom  size  for  R < 0.81.  This  modifi- 
cation provides  a physical  basis  for  the  composition  trends 
observed  in  ASDPs  for  the  most  stable  metallic  glasses.  Ap- 
plication of  this  modified  model  also  provides  a reasonable 
ranking  of  glass  forming  ability  in  complex  metallic  glasses 
[38].  However,  additional  work  is  still  required  to  allow  this 
model  to  rigorously  account  for  the  addition  of  several  so- 
lutes in  a single  alloy. 

Multicomponent  amorphous  alloys  that  contain  both 
larger  and  smaller  substitutional  solutes  relative  to  the 
solvent  are  generally  more  stable  than  binary  amorphous 
alloys.  This  feature  may  be  important  for  stabilization  of  the 
amorphous  structure.  Indeed,  as  substitutional  atoms,  the 
smaller  atoms  produce  compressive  lattice  strains  and  the 
larger  atoms  produce  tensile  strain  in  the  competing  crys- 
talline matrix.  These  opposite  strain  fields  may  attract  each 
other,  thereby  reducing  internal  stresses  and  forming  rela- 
tively stable  short-range  ordered  configurations  (clusters). 
If  these  clusters  can  be  systematically  linked  to  produce  a 
structure  with  long-range  crystalline  symmetry,  their  forma- 
tion will  stabilize  the  crystal  and  make  amorphization  more 
difficult.  However,  if  the  short-range  ordered  clusters  cannot 
be  linked  to  produce  crystalline  symmetry,  their  formation 
should  favour  amorphization.  Analysis  of  diffraction  data 
in  AI90Y10  and  A^YgNis  has  shown  that  there  is  an  ex- 
ceptionally strong  interaction  between  Y and  Ni  atoms,  and 
the  elastic  interaction  described  here  is  proposed  to  con- 
tribute significantly  to  this  behavior  [34].  The  model  based 


on  atomic  volume  difference  would  predict  a decrease  in 
the  glass  formability  of  alloys  with  both  larger  and  smaller 
solutes  relative  to  the  solvent,  and  so  contradicts  experience 
[37].  The  interpretation  provided  here  provides  a physically 
reasonable  explanation,  although  a more  rigorous  analysis 
is  required  to  validate  this  viewpoint. 

A completely  different  situation,  that  has  however  the 
same  nature,  can  be  considered  if  solutes  present  in  the  alloy 
are  small  enough  to  occupy  interstitial  sites  in  the  competing 
crystalline  lattice.  Indeed,  a higher  concentration  of  an  inter- 
stitial atom  is  required  to  reach  a critical  internal  strain  for 
destabilization  of  the  crystal  lattice  (amorphization)  when 
the  atomic  size  of  this  atom  decreases  relative  to  the  size  of 
the  solvent  for  R < 0.81,  as  discussed  above.  Producing  a 
tensile  lattice  strain,  these  interstitial  atoms  will  attract  sub- 
stitutional atoms  that  are  smaller  than  the  solvent  atom  and 
repulse  substitutional  atoms  that  are  larger  than  the  solvent 
atom.  In  the  former  case,  dense  and  stable  short-range  or- 
der configurations  may  be  produced,  which  may  stabilize 
the  amorphous  state,  while  the  latter  case  explains  why  bulk 
metallic  glasses  containing  large  amounts  of  small  interstitial 
elements  do  not  generally  contain  solutes  with  atomic  sizes 
larger  than  the  atomic  size  of  the  base  element  (solvent). 

A separate  approach  to  exploring  the  relationship  between 
the  concentration  of  atomic  species  and  the  structure  of 
metallic  glasses  has  recently  been  proposed  [19].  This  ap- 
proach is  based  on  the  concept  that  solute-centered  clusters 
may  form  common  structural  elements  in  metallic  glasses, 
and  that  these  elements  are  linked  together  to  form  the  glass 
structure.  Such  a concept  is  consistent  with  the  results  of  the 
previous  section,  where  specific  solute-to-solvent  radius  ra- 
tios favor  the  formation  of  atomic  clusters  with  a predictable 
value  of  N.  This  concept  is  also  consistent  with  the  stereo- 
chemically  defined  model  [1],  In  fact,  an  early  effort  to  con- 
struct a metallic  glass  structure  from  capped  trigonal  prisms 
showed  some  promise  [13].  In  the  approaches  described  here 
for  efficient  atomic  packing,  it  is  reasonable  to  expect  that 
a crude  relationship  may  exist  between  atomic  concentra- 
tion and  N.  If  a structure  is  comprised  of  solute-centered 
clusters,  and  if  each  solvent  atom  is  bonded  to  two  solutes, 
then  [19] 

CjMNqfei)  (3) 

where  C denotes  atomic  concentration  and  the  subscripts  i 
and  j represent  solute  and  solvent  atoms,  respectively.  This 
relation  is  satisfied  for  silicate  glasses  based  on  Si02  (C;  = 
0.33  and  Ny  = 4).  Further,  this  model  also  provides  a 
very  good  representation  of  metal-metalloid  glasses,  such  as 
Ni-B,  Fe-B  and  Co-P,  where  each  solvent  atom  is  coordi- 
nation with  ~2  metalloid  atoms  [39-41],  and  C;  ~ 0.18-0.2 
and  Nij  ~ 8.5-9.  In  Al-rare  earth  (RE)  glasses,  the  critical 
RE  atom  fraction  ranges  from  0.09  to  0.13  [42]  and  the  value 
of  Ny-ai  is  14.1  ± 1.5  [12].  Thus,  this  simple  model  gives 
the  solvent  concentration  as  0.63  < Cj  < 0.92,  which  spans 
the  experimentally  determined  range  of  0.87  < Cj  < 0.91 
[42],  Of  course,  as  N becomes  large  it  is  not  likely  that  each 
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solvent  atom  can  be  shared  by  more  than  one  solute,  so  that 
the  factor  of  1/2  in  the  right-hand  side  of  Eq.  (3)  may  re- 
quire modification.  In  fact,  the  value  of  A/ai-y  in  AI90Y 10  is 
1.6±0.2,  and  in  AlgyYnNis  is  1.2±0.2  [12].  Nevertheless, 
this  simple  model  is  a first  attempt  to  provide  a crude  ability 
to  relate  atomic  concentrations  to  the  structure  of  metallic 
glasses. 

5.  Concluding  remarks 

Structural  models  for  metallic  glasses  based  on  a dense 
random  packed  paradigm  are  unable  to  account  for  the  key 
experimental  observation  of  high  density  of  an  amorphous 
alloy  relative  to  the  density  of  the  same  alloy  in  the  crys- 
talline state.  A new  model  based  on  efficiently  packed  local 
atomic  configurations  has  been  proposed.  Several  key  fea- 
tures of  this  model  are  summarized  and  discussed  here.  A 
preference  for  specific  solute-to-solvent  atomic  radius  ra- 
tios predicted  by  this  model  has  been  validated  by  a strong 
correlation  with  radius  ratios  in  a wide  range  of  known 
metallic  glasses.  A dependence  of  solute  concentration  on 
atomic  size  is  also  discussed  in  the  framework  of  this  ef- 
ficient atomic  packing  model.  Finally,  specific  idealized 
atomic  configurations  are  being  constructed  that  are  con- 
sistent with  this  model,  and  validation  of  the  occurrence  of 
these  configurations  in  actual  structures  is  being  pursued  by 
consistency  with  available  experimental  and  computational 
results. 
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